Hybridization driven quantum critical behavior in weakly-itinerant ferromagnets 
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We investigate the unusual magnetic properties of nearly-critical, weakly-itinerant ferromagnets 
with general formula UTX, where T=Rh,Co and X=Ge,Si. As a unique feature about these systems, 
we show that changes in the Vdf hybridization control their proximity to a ferromagnetic instability, 
and determine the evolution of: the ground state magnetization. Mo, the Curie Temperature, Tc, 
the density of states at the Fermi level, N{Ef), the resistivity coefflcient. A, and the specific 
heat coefficient, 7. The universal aspect of our findings comes from the dependence on only two 
parameters: the Td bandwidth, Wd, and the distance between Td and U/ band centers, Ctj — Cuf. 



PACS numbers: 

In the past few years, a new class of materials has 
attracted significant attention due to their ability to sus- 
tain coexisting ferromagnetism and superconductivity at 
ambient pressure. These are ternary U compounds of 
general formula UTX, where T= Rh, Co, and X=Ge 
[TJ [2]. Together with UGe2, where superconductivity 
with Tc = 0.8 K is obtained under appUed pressure in 
a regime where ferromagnetism is still robust, Tc — 35 
K [5], and possibly Ulr [1], these materials have joined 
the special family of spin-triplet, weakly-itinerant ferro- 
magnetic superconductors. Furthermore, the observation 
that both superconductivity and ferromagnetism are very 
susceptible to applied pressure, indicates that the Cooper 
pairing is being induced by the proximity to a ferromag- 
netic quantum critical point (QCP) (5j. 

The search for a ferromagnetic QCP has thus encoun- 
tered a natural candidate, URhGe. Here the QCP could 
be reached through the chemical substitution of the tran- 
sition metal element, Rh, since other isostructural com- 
pounds are either paramagnets, as URuGe, or even closer 
to the ferromagnetic instability, as UCoGe. Thus it 
seemed natural to investigate the evolution of the fer- 
romagnetism and/or superconductivity in URhi_^Ta;Ge, 
for T=Ru or Co [5^. This was considered by Sakarya et 
al. in |6j, but, although the suppression of the ferromag- 
netism went according to expectations in the case of Ru, 
at Xc = 0.38, for the case of Co, instead, a huge and 
unexpected enhancement of Tq (larger than 100%) was 
observed in Tq, before it dropped down to less than 3 K, 
as a; — 1. The unexpected aspect of this finding comes 
from the fact that the isovalent doping in URhi-ajCoj^Ge 
does not introduce or remove carriers and the enhance- 
ment of Tc is much larger than what would be expected 
from volume compressibility effects (less than 30%) [SI. 

In this letter we describe the mechanism behind such 
unusual behavior of Tc as a function of isovalent disor- 
der, X, in URhi_3;Co^Ge, and we explain why a similar 



behavior is not observed in URhGci-ajSia;. We show that 
the evolution of Tc'(x) is governed by changes in the hy- 
bridization, Vdf, between the transition metal T=Rh,Co 
d-band and the actinide U /-level, which influences also 
other quantities such as density of states at the Fermi 
level, N{Ep), the T^ resistivity coefficient. A, and the 
speciflc heat coefficient, 7. As we shall demonstrate, hy- 
bridization ultimately controls the proximity to the fer- 
romagnetic instability, or QCP, where marked changes of 
behavior are observed for different ground state proper- 
ties in such weakly-itinerant ferromagnetic systems. 

The ternary URhGe and UCoGe compounds crystal- 
ize in the TiNiSi orthorhombic crystal structure with 
P{nma) space group. There are 4 chemical formula 
units per unit cell and each of the 4 U, Rh/Co, or Ge 
atoms occupy one of the 4(c) crystallographic positions 
(x, 1/4, z), (1/2 - X, 3/4, 1/2 + z), (1/2 + x, 1/4, 1/2 - z), 
and (— x, 3/4, — z), resulting in a primitive lattice. The 
relative positions of the U, Rh/Co and Ge atoms in the 
unit cell were determined by X-ray (300 K) and neutrons 
(20 K) [7 : xu = 0.9959(9), zu = 0.2038(05), XRh/co 
0.5716(7), ZRh/Co - 0.5730(17), xce = 0.6103(9), zce = 
0.5898(10), each of these quantities being measured in 
units of the lattice parameters a = 6.8552 A, 6 = 4.3274 
A, and c = 7.5010 A for URhGe, and a = 6.8450 A, 
b = 4.2060 A, and c = 7.2220 A for UCoGe. Samples of 
URhi_3;Co:,Ge with x = 0.0, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, and 0.9 were prepared by arc-melting the constituent 
elements (depleted U from Los Alamos National Labora- 
tory, 99.9 -I- % purity, Rh/Co/Ge from Johnson Matthey 
Aesar, 99.95%/99.995%/99.9999% purity) under a puri- 
fied Ar atmosphere. Samples were x-rayed using powder 
diffractometry in the University of Florida Major Anal- 
ysis Instrumentation Center. 

In a previous work on the effect of disorder on Tc [S\ , 
we found that sizable variations in the electronic mean 
free path, £, could drive the system from a clean towards 
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FIG. 1: (Color online) Experimental data for Tc as a function 
of doping in URhi_2;Coa;Ge (right vertical axis), as well as 
resistivity data (left vertical axis). 



a diffusive regime, where Tc was found to be enhanced. 
We thus wanted to compare the evolution of Tc (x) with 
the variation, with x, of disorder in URhi„2:Co3;Ge and 
URhGei_2.Sia;, using, as straightforward measures of the 
disorder, the residual resistivity pg or the inverse of the 
residual resistivity ratio, po/p(300K). As shown in Fig. 
[1] the Curie temperatures measured via magnetic sus- 
ceptibility for the as-melted chunks and the resistiv- 
ity bars (formed by sucking the molten material into a 
water-cooled mold, see [8J) agree well. However, as Fig. 
[1] also shows, the measures of disorder (either p(0) or 
p{0) / p{300K)) did not appear to vary as Tc{x), where 
the peak at a: = 0.6 is quite marked. The conclusion is 
that, for URhi_j;Coa;Ge, the electronic mean free path, 
£, is always large enough, kp£ ^ 1, such that the system 
is always in the clean limit, and thus the variation in Tc 
must be controlled by a different mechanism. 

The magnetic properties of ternary actinide com- 
pounds with general formula UT(Si,Ge), are determined 
by the SJ-nd-hybridization, Vdf, with n — 3,4,5 [9J. 
Quite generally, as one goes down a fixed column in 
the periodic table (isovalent elements), the 3d, Ad and 5d 
bands of the transition metals become broader, as well as 
lower in energy with respect to the position of the /-level 
of the U ion [9^. This is in fact observed for the UCoAl, 
URhAl, UlrAl series [9 , as well as for the UCoGe and 
URhGe compounds here under study |10l [TT] . As a con- 
sequence, the Vdf hybridization can be written as |12| 



df 



WdWf 



c. 



Uf 



(1) 



where Wd, Wf are the transition metal and /-level band- 
widths, and and Cu^ are the centers of the T^; and 
U / bands respectively. In terms of the muffin-tin orbital 
(MTO) theory by Andersen p2j, the above formula be- 
comes Vdf — (rjdfh'^ /m*)Jr^r'j:/d^, where m* is the con- 



duction band effective mass, rjdf are dimensionless an- 
gular momentum dependent coefficients, and r/ are 
the atomic radii of the nd and 5/ elements, and d is the 
interatomic distance f 14' . 

According to density functional theory (DFT) calcula- 
tions, for URhGe the 4(i-band of Rh starts at about —6 
eV, relative to the Fermi level, and extends beyond 2 eV 
above the Fermi level ^10) . thus being very broad. Fur- 
thermore, the Fermi level is located in a region of high 
DOS with 5/-character, see Fig. [2] The Stoner criterium 
is satisfied and URhGe exhibits ferromagnetic order with 
Tc ~ 9.5 K. This is also consistent with the high elec- 
tronic contribution to the specific heat for URhGe, which 
has been estimated as 7 « 115 mJ/mol (after sub- 

traction of the magnetic contribution). The results from 
DFT for URhGe also agree quite well with LSDA+U cal- 
culations [15j. For UCoGe, DFT shows that the Sd-band 
of Co starts at about —4.5 eV and extends up to approx- 
imately 1.6 eV, being, as such, much sharper |11] . Here, 
instead, the Fermi level is located in a region of much 
smaller 5/ DOS, see Fig. [2] and the Stoner criterium 
is barely satisfied. In this case, UCoGe is very close to 
the Stoner instability and is thus a nearly critical fer- 
romagnet with Tc « 3 K. Again this is consistent with 
the smaller electronic contribution to the specific heat for 
UCoGe, given by 7 w 65 mJ/mol [B]. 



URhGe 



UCoGe 












































































8 -U -4 


2 


2 4 



E - (eV) 



E - E^ (eV) 



FIG. 2: (Color online) Total DOS for URhGe and UCoGe. 
The blue (diamond) points are numerical data from Eqs. |4] 
and the solid black lines are numerical data from LSDA+U, 
for URhGe 153, and from DFT, for UCoGe [TT]- 

The microscopic mechanism for the unusual behavior 
for the Curie temperature in the URhi-^^Co^Ge series 
can now be described. As Rh Co, the d-bandwidth Wd 
becomes smaller as it changes its character from 4c? — 3d 
(we consider Wf fixed, for simplicity). The smaller Wd 
reduces the hybridizaiton, Vdf, which causes a large in- 
crease of the /-DOS at the Fermi level. Now, within the 
usual approach to itinerant ferromagnetism, one writes 
Tc (X (/7V(i;F)-l)^'=/''+^\ where J is the so called Stoner 
integral and z — 3 for a clean itinerant ferromagnet [16J. 
Thus, we see that an increase of N{Ep) leads to an en- 
hancement of Tc- On the other hand, as doping fur- 
ther increases, the bottom of the d-band is shifted to 
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FIG. 3: Mo from Eqs. Q. For increasing Vdf, a clear kink, 
separating the weakly and strongly hybridized regimes, occurs 
at around Vdf ~ 1.45, and is also observed experimentally |18| . 



higher energies, Ct^ moves towards Cu^ , since the 3c?- 
band is closer to the /-level than 4(i-band. The closeness 
between the d and / bands increases the hybridization, 
Vdf, which reduces the 5/-D0S at the Fermi level, reduc- 
ing Tc considerably. Thus, the non-monotonic behavior 
observed experimentally in the Tc{x) phase diagram of 
URhi_a;Coa;Ge results from the combination of the two 
effects on Vdf (or N{Ep)) when taken simultaneously. 

The stabilization of weakly-itinerant ferromagnetism 
in the periodic Anderson model has been studied in [171 



where t is the nearest neighbor hopping integral for 
the dispersing d-band, ^ are the position of the d, /- 
states, Vdf is the hibridization matrix element between 
the two, and U is the on-site Coulomb repulsion for the 



/-level, with n£. = flafic- For large enough U > Vdf 



a Hartree-Fock mean field decoupling can be justified. 



^ f f 
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where (n[^) gives the average occupation for the / level 
at the site i, including the two possible spin projections. 
We now introduce two Hubbard Stratonovich fields, Ni 
and Mi, associated, respectively, with the total number 
of / electrons, and the / electron magnetization, per site, 
= ('^{-^^ + ("4)' and Afi = {"^{^1 ~ {"'^{t) ■ We assume the 
system as homogeneous and replace Ni — rif — Nf/N, 
where Nf is the number of / electrons and N is the total 
number of sites. After this assumption we end up with 

H = Es^W<^^k,.+£/E/k'.^/k'.- + TE^^' 



Vdf 



]ifLdK. + - 1 E ^^^(4 - <) 

where ec(k) = 2i J2i=i {kia) is the dispersion, and 
the renormalized /-level is shifted toef—ef + Unf/2. 

We introduce ferromagnetic order by writing Mj — Mq. 
The Hamiltonian becomes H = NUM§/A + Y,\^C'^T-iC, 
where we have defined C"^ — (c|^ ^cj^. ^/j^ ^/j[ ^) and 



k 



n = 



e,(k) Vdf 

ec(k) Vdf 

Vdf Ef-h 

Vdf ef + h 



(3) 



The above Hamiltonian is quadratic in the fermion opera- 
tors and can be easily diagonalized. The new eigenvalues 
are (/ = lower ,u = upper branches) 

Eu,±{y) = \ [(ec(k) +ej±h) + y^(e,(k) -Sf^ + 41/|J 
£,,±(k) = 1 [(e,(k) + Ef ±h)- y^(e,(k) -SfT + 41/|J 

where h — UMo/2. The self-consistency equations for 



Mo and rif are given by Mq = (l/N) J^kii'^U^ " ('^k;)} 
and rif — {1/N) X^kK^-kt) + ("'H'' J' '^^ explicitly as 



= ^EE^^g^^(^) 



k <T=± 



nFD{Eu,aW - /i) + — TTT 77-r nFD(El^a[ii) ~ m) 



TV ^ ^ 



k a=± 



i;„,,(k) (k; 

£;„,^(k)-£,(k) 
i;„,,(k) (k 



£;,,,(k)-s„,,(k) 



(4) 



where nFD{E - fi) ^ (e^(E~t^) + ^jth /? = l/fc^T. 
We find Mq{x = 0) = 0.22^B/U-atom, for URhGe, and 
Mo{x = 1) = 0.06/XB/U-atom, for UCoGe, in very good 
agreement with experiments [18 , showing that in UCoGe 



Vdf is indeed much larger, see Fig. |3] We have also cal- 
culated the total density of states for the above band 
structure using Eqs. 4] with n/ = 0.4 [1] and J7 = 9 (in 
units of t). The results in Fig. |2]are in very good agree- 



4 



ment with LSDA+U HU and DFT data [TT]- Notice that 
the spHtting between + and — states near Ep is propor- 
tional to h = [/Mo/2, being thus much larger for URhGe 
than for UCoGe, see Fig. |2] The change in behavior of 
M^oiVdf) occurs exactly when Vdf becomes larger than 
such splitting, and Mq decreases linearly as a function of 
Vdf, before it vanishes for V^^ « 1.9 (in units of t), at the 
ferromagnetic QCP. 

Since the d-bandwidth, Wd, decreases monotonically as 
Rh is replaced by Co, we parametrize 



Wdix) 



Wd^'^{l-x) + W^''x, 



(5) 



where the values for W^'^ and agree with pUl [TT] . 

giving t^^ = 1 eV and — 0.6 eV. Analogously, 
= Cxa approaches continuously to e° = C'cz/i with 
C^^ = -2 eV, C^^ = -1.5 eV, and Cuf = eV, before 
hybridization. However, this approach, though mono- 
tonic, deviates from linearity [9 and thus we write 



ACdfix) = AC§''{l-x)+AC^/x+6'x\l~x)+S"x{l~xf, 

(6) 

where d' and 5" are adjustable parameters. Below we plot 
Tc{x) for URhi-a^Co^^Ge where the Vdf hybridization in- 
duced, non-monotonic behavior for Tq can be observed. 
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FIG. 4: Plot of Tc{x) from Eqs. Q and of A{x) from Eq. 
([7|. Experimental data for A was extracted from |18) . We 
have fixed Wf so that Tc{x = 0) = 9.5 K, in which case 
Vdf{x = 0) = 1.42 and Vdf{x = 1) = 1.73, see also Fig.|3] The 
values for the adjustable parameters, S' = 1 and S" = 0.8, 
where used throughout the paper (in units of t). Notice the 
kink in Tc{x), for x — 0.9, at the Vdf crossover, see Fig. [s] 

The nontrivial evolution of N{Ep) governed by Vdf 
has also measurable consequences in transport proper- 
ties. For example, the resistivity of URhi_a;Coa;Ge be- 
haves as p ~ AT^, for all x, where A, the resistivity 
coefficient, exhibits also a nonmonotonic behavior as a 
function of x jTB]. Although it is generally difficult to 
separate electronic and magnetic contributions to A, for 
the case of weakly-itinerant ferromagnets, the dominant 
magnetic scattering contribution to the resistivity is [16| 



Pmag ^ \INiEF) - l\-^I^T'' = AT 



(7) 



In fact, A can be thought of as a direct measure of 
the proximity to the ferromagnetic QCP, being large for 



URhGe but even larger for UCoGe, smaller at interme- 
diate doping < a; < 1 |TH|. Thus pmag should have the 
opposite behavior as Tc, see Fig. |4] right panel. 

It is interesting to compare the above results for 
URhi_2,Coj;Ge with the case of another isovalent- 
disordered compound, URhGei_a;Sia;. The substitution 
Ge — > Si in URhGe i_j:Sia; does not modify the Vdf hy- 
bridization, and causes, at most, a small broadening of 
the Rh, 4c?-band, due to substitutional disorder. Such 
small broadening causes an almost negligible reduction of 
the 5/-D0S at the Fermi level, and thus Tc is expected 
to be slightly reduced at the point of largest disorder, x — 
0.5. This is exactly what we measured for URhGeg.sSio.s 
where we found Tc{x = 0.5) = 8.5 K, only 10% smaller 
than the endpoints Tc{x 0) = Tdx = 1) = 9.5 K. 

We have shown that in nearly critical, weakly itin- 
erant ferromagnets, such as URhi_a;Co2:Ge, the prox- 
imity to a ferromagnetic QCP is controlled by the hy- 
bridization, Vdf, between the transition metal and the 
actinide. Smaller values for Vdf drive the system away 
from the QCP, while larger values for Vdf pushes the sys- 
tem towards the QCP. The unusual behavior observed 
for Tc{x), Mq{x), and A{x), can all be quantitatively 
understood from the nontrivial evolution of Vdf (x) as the 
4c?-band of Rh changes character to the Sd-band of Co. 
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